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Abstract 

The bound-state problem for the pion as a quarkonium with the fun- 
nel (Coulomb-plus-linear) interaction is solved in a framework that com- 
bines the bilocal approach to mesons with the covariant generalization of the 
instantaneous-potential model. The potential interaction leads to dynam- 
ical breaking of chiral symmetry. However, the Coulomb potential leads to 
ultraviolet divergences that must be subtracted. A careful choice of the renor- 
malization prescription is needed in order to get the correct chiral limit. The 
mass, the lepton decay constant of the pion, as well as the pion decay width 
in two photons are calculated. 
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I. INTRODUCTION 



The activity in the field of relativistic bound states has recently increased. References 
[ t]-[5|] are just some of examples. These papers were motivated, among other things, by the 
desire to formulate a covariant treatment for quarks interacting through a potential which 
would hopefully mimic nonperturbative QCD. 

Our work is a continuation of the line of research where interquark interactions were 
modeled by the instantaneous potential. For the case of light quarks, Refs. may serve 
as paradigmatic examples of such calculations. A related approach, using the Nambu-Jona- 
Lasinio model, can be exemplified by Refs. |§. 

The successes of the potential model in describing heavy quarkonia are well known. For 
light quarks, however, besides many important qualitative successes, the potential model 
also exhibited weaknesses: first, the uncertainty concerning the question what potential 
can provide a realistic and yet reasonably tractable interaction without many free parame- 
ters; and second, the noncovariance of the instantaneous-potential approach to such highly 
relativistic constituents as u, d, and s quarks. 

As an example on the successful side, Le Yaouanc et al. |J demonstrated the appear- 
ance of dynamical (spontaneous) chiral-symmetry breaking by generating the dynamical 
quark mass as well as the pion as a Goldstone boson in the chiral limit. They studied the 
power-law interaction (V oc r a ) between massless quarks and antiquarks most exhaustively 
in the simplest case of harmonic potential (V oc r 2 ) where the gap equation (Schwinger- 
Dyson equations) reduced from an integral equation to a differential one. Adler and Davis 
0] formulated the renormalization procedure for the Coulomb-like potential and performed 
a concrete numerical calculation for the linear confining quark-antiquark potential. In the 
latter case, only infrared divergences appeared. Trzupek || applied their renormalization 
procedure to the more realistic case of the funnel (linear-plus-Coulomb) potential where 
ultraviolet (UV) divergences were present. This situation was further complicated if finite 



quark masses were present [ 1 C . 1 1 ] . Nevertheless, the aforementioned weaknesses led to un- 
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satisfactory quantitative results. For example, the investigations of Refs. g-§ could not 
yield the value for the pion decay constant better than four to five times smaller than the 
experimental one, when the meson spectrum was fitted correctly. This was assumed to be the 
consequence of noncovariance |J. (Furthermore, it was expected that finite current quark 
masses would improve the results). Our attention was thus attracted by the covariant gener- 
alization of the instantaneous-potential approach to bound states, formulated by Pervushin 
and collaborators []12-15|. This approach was first applied to the harmonic potential as the 



simplest case ||16|| : however the first concrete and correct numerical results were obtained 



by our group |1y|| . Besides covariance, the effect of the finite current quark masses was also 



included |T7| , while previous investigations in this line of research 0-§J had been concerned 



with massless quarks only. Therefore, they had not been pertinent for studying the depen- 
dence of the pion mass on the model parameters since in the chiral limit the pion mass is 



vanishing. The mass of the pion from our Ref. |17j behaves as the square root of the current 
quark mass, which is precisely the correct behavior of the (pseudo) Goldstone boson (see, 
e.g., [|TBJ). The pion decay constant, however, was found in Ref. [[T7]] to be F n « 35 MeV for 
the harmonic-potential strength which reproduced the experimental pion mass. The result 
turned out to be practically the same as in Ref. §, F n = y/2f„ = ^2 20 MeV « 28 MeV. 
The covariant approach removed certain ambiguities present in the definition of the pion 
decay constant in Ref. 0, but did not solve the problem of its too small value. Obviously, 
further studies of the form of the interaction potential were needed. In the present work we 
therefore examine the funnel (Coulomb-plus-linear) potential 

4 a 

V(r) = V c {r) + V L (r) = -(-- ^ + ar) , (1.1) 

3 r 

since it is known that, in QCD, the short-distance interactions are dominated by the 
Coulomb interaction, while in the long-distance (small momentum, or k — > 0) regime, 
a s times gluon propagator seems to behave as 1/fc 4 , corresponding to a linear confining 
potential in the coordinate language. However, the ultraviolet divergences caused by the 
Coulomb part pose some new difficulties. In fact, the issue of renormalization of the bound- 



state equations for quarkonium with instantaneous interaction is the main point of this 
work. After sketching in Sec. [TI| how the representation of mesons by bilocal fields leads 
to the Schwinger-Dyson equation (SDE) and the Bethe-Salpeter equation (BSE) in ladder 
approximation, in Sec. |TTT| we formulate a renormalization scheme for the SDE with the 
funnel potential. We discuss the limitations which such a scheme must suffer when various 
approximations are introduced. We also compare our renormalization procedure with the 



ones used so far in this context. In Sec. [TV] the Salpeter equation for the pion is solved, and 
in Sec. the pion decay constant is obtained. In Sec. [Vl| we calculate the 7r° — > 77 decay 



width and conclude in Sec. VI]. 



II. MESONS AS BILOCAL FIELDS 

When trying to model nonperturbative QCD, one may consider a very general interaction 
kernel K entering in the effective action: 

+ 2/ d Sj^aiy)qpix)[K{x-y)\ auPl . a2A qply)q a lx)} . (2.1) 

where m is the current quark mass matrix, m = diag(m u , m<j, m s ). oti and /3j are spinor 
indices, whereas color indices and flavor indices are suppressed. In fl2.1|) the summation over 
repeated indices is understood. We assume that the interaction kernel K(x — y) can lead to 
a bound qq system. In Q2.1J) we have introduced L(x), an external operator coupled to the 
quark current. For example, it can be the leptonic current l^{x): 

L(x) = ^W 1 -^ , (2-2) 

or a photon A^x): 

L(x) = eA>*(x)^ = e 4(x) . (2.3) 

Such external operators will make possible the weak and radiative decays, but the internal 
structure of hadronic bound states will be dictated by the model kernel K. 
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One can construct a theory of meson bound states by eliminating bilinear structures 
Qa(x)Qp(y) in favor of bilocal fields Xa,p{ x iV) IHhH' (introduced through the path integral 
in the generating functional) and then integrating out the remaining quark fields. In this 
way the action ( |2.1| ) becomes pi 



W eff [q,q\ -> W eff [x\ 

= iN c Tr\n[ift-m — L — x] + -^(Xi^x) 

= N c {iTrln(ip-m) - tTr £ -{{ip- m)~\L + + ~(x, K ~ l X)} , (2-4) 

n=l 71 2 

where we have suppressed all indices and used shorthand: 

{.X,K~ l x) = J d 4 xd A yXf3 ia2 (x,y)K-^ i . a2:f32 (x,y)x l 3 2 a 1 (y,x) , (2.5) 

and where Tr (with the capital "T" ) also includes the integration. (Below, "tr" , with small 
"t", will denote a trace not including integration.) We can drop the external (weak or 
electromagnetic) operator L while studying the bound-state equations determining the in- 
ternal hadron structure. We shall reinstate L later, while studying weak and electromagnetic 
decays. 

We determine the classical solution xo conveniently written as Xo( x ~ y) 
= XX^) v) ~ i : n5 A {x — y), by varying W e ff with respect to x : 

6W f [x] = . (2.6) 
OX 

This yields the SDE for the quark self-mass operator E in the ladder approximation: 

S(x — y) = m5 4 (x — y) + iK(x — y)S(x — y) , (2.7) 

where the "dressed" quark propagator is defined by 

S^(x — y) = (ip — m)8^ 4 '(x — y) — Xo 

= ipS {4) (x-y) -E(x-y) . (2.8) 
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Next, we expand the fields x i n the action around the minimum, x( x >y) = Xo( x jV) + 
Ai(x, y) = J2{ x ~ y) ~ ^^(x — y) + Ai(x — y). As it will turn out that the fields A4(x, y) 
represent mesons, we separate the part of the action containing A4(x, y): 

W ef f \x\ = W eff [xo + M} = W eff Ixo] + Weff [M] , (2.9) 

where this part is given by 

~ N 00 1 

Weff[M] = -f(M, K- l M) - iN c J2 -Tr$ n (2.10) 

Z n=2 n 

and where we have introduced 

$(x,y) = Jd 4 zS(x,z)M(z,y) (2.11) 

and, as before, Tr$ n denotes 

Tr$ n = tr J d i xid A x 2 ---d' >t x n ^(xi,X2)^{x2-,x^)...^{x n ,Xi) . (2.12) 

Up to the terms of the order 0(A4 2 ), varying the action W e //[xo + -M] or ; equivalently, 
We//[A^], with respect to the fluctuations A4(x,y) gives the BSE in the ladder approxima- 
tion for the bound state of a quark and an antiquark whose spectra and propagators are 
determined by Eq. (|2.7|) and whose flavors a, b are now explicitly written out: 



M ab (x, y) = iK(x - y) / dx'dy'S a (x - x')M ab (x', y')S b (y' - y) . (2.13) 



Equation ( 2.13|) is the BSE in the somewhat improved ladder approximation, since the quark 



propagators in it are not the free, bare ones, but S, containing the nontrivial self-energy 
function E. 

In this paper we work in the isosymmetric limit, being concerned with pions only. Not 
having to distinguish quark masses for different flavors allows us to simplify the notation 
in the rest of this paper by dropping the indices a and b. (However, in processes involving 
kaons, it will be necessary to keep track of quark flavors carefully.) 
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III. CO VARIANT GENERALIZATION 
OF THE POTENTIAL APPROACH AND THE RENORMALIZED 
SCHWINGER DYSON EQUATION FOR THE FUNNEL POTENTIAL 



To solve the SDE and the BSE in practice, one must limit oneself to a tractable interaction 
kernel K, which is given by 

K(k) = -tC F g 2 Y ® YD^(k) . (3.1) 

where Cp is the second Casimir of the quark representation, here 4/3 for the case of SU(3) 

triplet, g is the strong coupling constant, and D is the gluon propagator. An instantaneous 

approximation to the kernel K leads to the potential model, 

~ ~ Wk 1 

K(k) « z 7 ° ® -y°V(k) - ij j ® j l V T (k) [5 jl - —} . (3.2) 

\k\ 2 

Further approximation consists in neglecting the transverse gluon exchange. Thus, the 
kernel K becomes 

K{k) w ij° (g) 7°V"(fe) . (3.3) 



Pervushin and his group f2J[I^-|T5| have found that the covariant generalization of the 
potential approach is possible, provided that the kernel K is of a special form K n : 

K n {k) =ifi®fiV(k ± ) , (3.4) 

where n is the timelike unit vector in the direction of the total momentum of the bound sys- 
tem, = P M / \f~P^. For any vector k^ 1 = k^ + kj_, the components parallel and perpendicular 
to this axis are 

E = n^kp, k P = k-n = k- P/VP*, 



(3.5) 

jfcf = k» - ktf, k i ■ P —— . 



The relativistic covariant formulation of potential models, given by the kernel 
guarantees that the correct dispersion relation for the momentum and mass of the bound 



state, P 2 = M 2 , is fulfilled. 



For the interaction kernel of the form ([3.4Q , the SDE (|2.7|) Fourier-transformed to mo- 
mentum space is 

S- 1 (p)=t-™ + ij ^LfS(kW(p ± ~ fc±) • (3-6) 

Equation ( |3.6| ) is valid for an arbitrary reference frame. For the quarkonium rest frame, 
n = (1, 0, 0, 0), and in the chiral limit, fl3.6|) reduces to the SDE studied by, e.g., Le Yaouanc 
et al. [f| in their noncovariant approach with the power-law interaction, V(r) oc r 7 . In 
Ref. [17| we have already studied Eq. (|3.6| ) for the particularly simple, harmonic case, with 



7 = 2, where the integral SDE reduces to differential equations. 

The UV divergences due to the Coulomb part of the potential require renormalization 
and introduction of counterterms. This will change Eq. ( p.6|) and its rest-frame version. 
Following Ref. JF| and its generalization to the massive case [|IT)| , |lT|| , we use the equations 



for renormalized vector, axial-vector and pseudoscalar vertices, and Ward identities, to set 
the renormalized SDE in the ladder approximation, 

S-\p) = Z4 - Z m m - ig 2 C F J ^-^s^D^p - k) , (3.7) 

where Z 2 and Z m are the wave function and mass renormalization constants defined by 

S = Z 2 S and m = % CL m , (3.8) 

^2 

where So and m are the bare quark propagator and the bare quark mass, respectively. 
Neglecting the retardation effects, Eq. (|3.2|), and the transverse gluon exchange, Eq. (|3.3|), 
in (|3.7|) we are provided with the renormalized version of Eq. ( |3.6| ): 

S-\p) = Z 2 i> -Z m m + ij -0-^Sik^Vip - k) . (3.9) 

Since the perpendicular part of a four-vector reduces to the corresponding three-vector 
in the quarkonium rest frame, e.g., k± — > k, we will use the noncovariant notation to the 
end of this section. This will make easier the comparison of our results with the results of 
other authors in this line of research. Of course, all the expressions can be generalized back 
to those valid in moving frames, by substituting k —>■ k±, k° — > kp, and 7 — > fi 
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Let us demonstrate the multiplicative renormalizability (MR) of this equation. The 
renormalization of the product g 2 D is 

dD = AyD, (3.10) 

^2 

where the subscript refers to the bare quantities and Z\ is the vertex renormalization 
constant. The gauge invariance implies Z\ = Z 2 and the renormalization-group (RG) in- 
variance of g 2 D. However, the ladder approximation is consistent with Z\ = 1, i.e., with 
no vertex renormalization (see, e.g., fl2"5|). So, the renormalization of g 2 D, and hence of V 



should be 

9 2 D = (U 2 g 2 D , (3.11) 

^2 

V = (^) 2 V. (3.12) 

^2 

Now, suppose that {Z 2 , Z m } and {Z' 2 , Z' m } are two sets of renormalization constants. 
They may correspond to two different renormalization scales fi and //. Since there is a 
definite relationship between the bare and renormalized quantities, we know the relationship 
between the quantities renormalized by primed and unprimed Z's. Concretely, using ( p.8|) 
and ( [3.12j ), we transform the SDE, Eq. ( p.9|) , and find that it does not change its form, so 
that the MR holds. Before further discussion of the MR, we shall rewrite the SD equation 



( p.7| ) using the conventional ansatz of the quark propagator S through the functions ui and 

(p: 

S~\k) = A;°7° - oo(k)C 2 (k) , (3.13) 
where the matrix ( is defined as 

C(fc) = sin ^y2(fc) — k ■ 7 cos ^<p{k) . (3.14) 
We can express the quark propagator as 

-(1+1°) ifl-O' ) 

S{k) = -C(fc)[ /7 2 ,„ . + „ 2 / n } . ]C(fc) • (3.15) 
v ' y n u(k) - k° -is u(k) + k° -te 1 
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Inserting ( |3.13|) into the renormalized SD equation ( |3.7| ) yields the following integral equa- 



tions for uj and <p: 

1 f d?k 

uj(p) sin (p{p) - Z m m + -J J^A3 sin ^ k ^ V ^ p ~ fc ) = , (3.16a) 

1 f d 3 k - ~ 

uj(p)costp(p) - Z 2 \p\ + -J -— ^ cos (p(k)(k-p)V(p -k) = 0. (3.16b) 

Additionally, the equation (Z 2 — l)p° = arises, i.e., Z 2 = 1. However, the Eq. ( |3.16b| ) 



demands that Z 2 to be an infinite constant. Adler and Davis J7| have resolved that contra- 
diction by splitting Z 2 into two parts, 

Zrf^ Z oP V - Zp ■ 7 , (3.17) 

and setting Zq = 1. However, we have to remind ourselves that Z 2 is the wave-function 
renormalization constant, which defines the renormalization of the quark propagator, So = 
Z 2 S. We are faced with a dilemma, namely, whether S should be renormalized with Zq 
or with Z. As one can expect, both possibilities change the form of the renormalized SD 
equation Q3.9p. Obviously, the consistency of MR is violated because one has to split Z 2 into 
Z Q and Z. If unprimed and primed renormalization constants correspond to two different 
renormalization scales \i and //, respectively, this inconsistency automatically shows that the 
invariance with respect to the changes of the renormalization scale \l is lost, and we cannot 
use the renormalization group (RG) equations to relate results for one arbitrary scale fi to 
that for some other scale fi' . Nevertheless, as noted by Brown and Dorey |[25|| , who explored 



the consistency of the MR of the SDE when various approximations are made, this does not 
mean that treatments that do fail such a consistency test cannot be useful, merely that is 
more difficult to relate their solutions to real physics. Of course, we must investigate the 
scale dependence of our results. We shall return to this point later in the text. 

A convenient renormalization prescription that determines Z and Z m uniquely is given, 



e.g., in Refs. [S6|,27|. The authors specify that the quark propagator S(p), for a given 
spacelike p 2 = —fi 2 , agrees with free-field theory. We adopt this choice, adjusted for a 
special form of the kernel (|3.4j) and the propagator (|3.13|) , 
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S-\p)\\ p \ = fl =^-m. (3.18) 



Imposing (|3.18|) on the SD equation ( |3.7p yields the renormalization constants 



1 r d 3 k - ~ 

Z = l + -jj^- 3 cos V (k)^-k)V( f i-k) , (3.19a) 

1 r d 3 k ~ 
Z m = l + —j sm V {k)V{p,-k) . (3.19b) 



Using (|3.19|) , the SD equation (|3.16| ) becomes 



u>(p) sin <p(p) — m 

1 r d 3 k 



r (I k ~ ~ 

+ 2 J ^ [ViP ~ ^ ~ Vifl ~ ^ = ° ' (3 ' 2 ° a) 

uj{p) cos (f(p) — \p\ 

+ ~ / -0-, cos ip(k)[(p ■ k)V(p - k) - ■ k)V(iM - k)] = . (3.20b) 



For m = 0, there is no mass renormalization, and the functions tpz> and ujd of the 
dynamical quark propagator Sp have to satisfy the equations 

1 r d 3 k 

u D (p) sin Lp D (p) + - J j^-^sm(p D (k)V(p- k) = 0, (3.21a) 
u D (p) cosy? D (p) - \p\ 

+ 2 y (2^p cos^(fc)[(p • fc)V(p - k) - ■ k)V{n - k)} = 0. (3.21b) 

Now, it is obvious that the SD equation (|3.20|) does not reduce to Eq. (|3.21|) in the limit 
m — ► 0, because lim m ^ Z m m ^ 0. Moreover, we see that <p(p) = <Pd(p) and ui(p) = ujd(p) 
is a solution to the SD equation ( |3.20|) for m = w! = u D (n) sm.(p D (p). So, the chiral 



limit is reached in the SD equation ( |3.20|) for m — > m' and not for m — > as it should be. 



This is an artifact of the renormalization prescription Q3.18] ). Pagels |2£| showed that the 
normalization condition (|3.18| ) precluded the presence of a dynamically generated term in 
the quark propagator S(p), and he argued in favor of Weinberg's zero-mass renormalization 
scheme. However, it is possible to recover the proper chiral-limit behavior by redefinition 
of the mass renormalization constant, Eq. ( |3.19b| ): 
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1 1 r d 3 k ~ 

Z m = l + —u D (fj,)sm(p D (fj,) + —j -^-^sm(p(k)V(fi - k) . (3.22) 



This corresponds to the redefinition of the normalization condition (|3.18|) : 



S 1 (ti\\p\ =/i= $- m -uv(p)sm(p D (p)]\ p \ =fi . (3.23) 
Equation ( |3.20a| ) becomes 

uj(p) sin (p{p) — u)d{i*) sin V9 D (^) — m 

+ \Swf ^ v{k)[V(p ~ k) - v ^ - k)] = ' (3 - 24) 

whereas Eq. ( |3.20b| ) remains unchanged. 



The renormalization constants Z, Eq. ( |3.19a| ), and Z m , Eq. ( |3.22| ), are defined in terms 



of the potential V. We are considering the case of funnel potential, which is a sum of the 
Coulomb-like potential Vc and the linear potential Vj> The integrals involving the linear 
potential are finite, so we can drop Vl from the definition of the renormalization constants, 
i.e., we can define 

1 r d 3 k - ~ 

Z = 1 + 2UT\J J^f cos ^ k ^ • - k) , (3.25a) 

1 f d^ k I f d^k 

Z ™ = 1 - ST 7 ^i-^W^ - *) + ^sinmVafp - k) . (3.25b) 

The first integral in the Z m definition, Eq. ( |3.25a| ), multiplied by m, now plays the role of 
m! . Omitting this term will cause the incorrect chiral limit behavior we have met using the 
renormalization constants ( 13.19 ). The renormalization constants ( |3.25| ) are those we have 



actually used in our numerical calculation. The corresponding SDE is 

1 f d 3 k 

uj(p)smp(p) - [m - - J j—^Wi(p D (k)Va{p-k)] 

1 f d 3 k ~ 
+ 2J J^f sin ^( fc )[^(P " k ) ~ V c(» ~ k )\ 

1 f d 3 k ~ 
+ 2~J J^ s ^^ k )Mp-k) = 0, (3.26a) 

u(p) cosy?(p) — \p\ 
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1 r d 3 k ~ ~ 

+ 2 J (2^p C ° S ^ (fc)( ^ ' k)VL ^ P ~ k) = ° ' (3 ' 26b) 



Let us now make a comparison with the renormalization prescription of Refs. P,|10|,|29 
The renormalization constant Z used by them is 

1 r d 3 k - ~ 

z = 1 + 2wf(M^- k)Vcip - kh (3 - 27) 

One way in which Eq. ( {3.27 ) differs from our Z, Eq. ( p.25b|) , is the absence of the factor 
cos </?(&). This is not a critical difference, as the modification of Eq. ( |3.25|) by making 
the substitution cos<p(k) — > 1 would change Z only by a finite value. As noted by Adler 
and Davis 0, the crucial point is that Eq. ( |3.27| ) corresponded to a momentum-dependent 
subtraction of UV divergences when the coupling constant a s was running, i.e., when a s 
was momentum dependent. (For that reason, they rejected the SDE with the Coulomb- 
like interaction of Finger-Mandula ||30|| .) For a s = const, however, they noted that the 
infinite part of ( |3.27| ) defined a momentum-independent UV subtraction. Nevertheless, this 
counterterm still has a momentum-dependent finite part. This is conceptually objectionable 
even if the UV infinity is successfully subtracted. 

The other possibility considered by Adler and Davis for the Coulomb-like countert- 
erm, again in connection with the unrenormalized equations of Amer et al. |3ll , |32| . This 
alternative Z is given by Eq. (2.17) of Adler and Davis and is obviously a completely 
momentum-independent choice. It can be shown that it is given by (|3.27|) when \p\ tends 
to zero: 

1 r rfik ~ ~ 

z = 1 + ^\Iwr (t " k)Vc{p - kh (3 - 28) 

It is evident that ( |3.28| ) is momentum independent simply because |p| is fixed to a specific 

value - zero. However, there is a subtlety that seems to have been unnoticed so far. After 

performing angular integration in Eq. (|3.28 ), we obtain 

Aa s f°° dk , , 

Z = 1 - — - / — . (3.29) 



3tt 
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This means that, for A; — > 0, the UV renormalization constant Q3.28| ) introduces a new IR 
divergence, which has not been present in the original SDE! If a s is not constant but runs 
with k, then a s have to be under the integral; however, this does not change the situation. 
One possibility for avoiding troubles for p — > is to introduce an IR cutoff, but this requires 
a new unphysical parameter. Alkofer and Amundsen [3~3~ used the renormalization constant 
( [3.281 ) but no IR cutoff, which makes their results suspicious Q. On the other hand, our 
renormalization constant (|3.25| ) also reduces to (|3.28| ) as — > apart from the factor 
cos ip(k) in the integrand. For \k\ — > 0, tp(k) — > n/2 and cos (p(k) — > 0, which makes the 
integration in ( |3.25|) IR regular. 

To summarize, our Z differs from that used (so far) most frequently in this line of research 
|7|JT0|j8]j29[] , in that the variable momentum is replaced by the fixed renormalization scale fi. 
At the same time it encompasses the alternative possibility [|7|,|3~3f for Z in the limit \x — > 0. 
The finite part is chosen so that no new IR divergence arises in this limit. 

Our solutions are given in Fig. [I], while numerical methods are detailed in the Appendix. 

The linear potential Vi makes the integrals in the SDE ( |3.2(j| ) IR divergent. The formulae 
( |3.2(j| ) can be rewritten such that the equations for tp and u are separated. The equation for 
ip needs no IR regularization The other equation expresses w as a functional of ip and 
requires IR regularization. We adopt the procedure used by Le Yaouanc et al. f|, where 



4 4,2 e-« r - 1 + £r 
—ar = -a lim — 

3 3 ?-o£2 



V L {r) = -ar = -a lim ^ , (3.30) 



which gives the IR-finite quark energy u(k). For light quarks, u(k) is negative for small 
momenta, as expected from earlier works, e.g., P]T6|j. This is in fact a signature of dynamical 
chiral-symmetry breaking. 

To conclude this section, we have successfully solved the SDE for massive quarks with the 
funnel-potential interaction, treating IR and UV divergences carefully. With the solution for 



1 Alkofer and Amundsen studied the temperature- and momentum-dependent a s , but the claim 
of the singularity in Z for p — > still holds. 
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(pip), we can proceed to solving and examining the Salpeter equation for a light pseudoscalar 
quark-ant iquark bound state. 

IV. THE SALPETER EQUATION FOR THE PION 

Solving the SDE ( |3.26| ) for <p{k) provided us with the necessary input for solving the 
BSE (|2.13|) . By Fourier transforming to momentum space, this equation becomes 

r(g|P) = i J 0^K(q - q')S(q' + |)r(g'|P)5( 9 ' - |) , (4.1) 

where 2q = p a — pb and P = p a + Pb are the relative and the total momentum of the quark- 
antiquark pair, respectively, and T{q\P) is the momentum vertex function of their bound 
state. Using the instantaneous potential for the interaction kernel reduces the BSE to the 
Salpeter equation. This equation is still manifestly Lorentz covariant for the special form 
( P-4| ) of the instantaneous interaction, as shown by Refs. [|1~3|-|T5|| . It is convenient to write 
this equation in terms of the quarkonium wave function ^{q±) or its transformed mate ip p : 

* P (<Z±) = i J j^c[S(q + j)T(q ± \P)S(q - £)] (4.2) 

In this work we are interested only in pseudoscalar mesons. It can be shown that Pip p = 
—if) p P, so that the decomposition of ip p for pseudoscalar mesons in the Dirac matrices is 
simply 

?P P (q±) = 75[£i(<Z±) + -7^H<1±)] ■ (4-3) 

v P 

The BSE ( |4.1|) written in terms of (^4.2|) and ( |4.3|) and for the form ( p.4|) is boost invariant. 
To solve it for L\ and L 2 , we are free to choose the rest frame of the bound system q^ = (0, q), 
P = [M n , 0). The mass M n of the pseudoscalar bound state, the pion, is the eigenvalue of 
the BSE. (The treatment of the kaon is in principle identical to that of the pion, except that 
one of the quarks would have to be the significantly heavier strange quark). 
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Inserting ( |4.2| ), ( f4.3|) , and ( |3.4j) in ( |4.1|) yields the coupled system of integral equations, 
which at first sight seems to be IR infinite, 

M^L 2 (p ± ) + 2u( P± )L 1 ( P± ) + J j^V( P± - k ± )In{k ± ) = 0, (4.4a) 

V(p± - /c±)[sinv?(p±) smip(k±) 
- (p± ■ k±) cos ip(p±) cos (p(k±)]L 2 (k±) = 0. (4.4b) 

However, as in the SDE, the IR infinity of the unregularized quark energy u(k) cancels 
the IR infinity of the integrals with kernels (see |?J). 

The numerics for solving (|4.4j ) is discussed in the Appendix along with the numerics for 
the SDE. The solutions Li and L 2 are displayed in Fig. [[]. 

Having solved the massless version of the SDE, Eq. (|3.26| ), we have found that 
ujd{h) sin </?£>(//) — > as /i — > oo. So, in the limit of infinitely large renormalization point 
/i, the normalization condition ( |3.18| ) approaches the normalization condition ( |3.23|) . The 



improper behavior of the theory, artificially generated by the normalization condition fl3.18| ) 



will disappear for — > oo. Fig. shows the pion mass as a function of the renormal- 
ized quark mass m for three different renormalization prescriptions. The dotted lines relate 
to the normalization condition ( |3.18| ), slightly modified by leaving out the linear potential 



from definition of the renormalization constants, and fi = 1 GeV. In this case, M v — > for 
m -> m'{fi) = u D {n) sin^ D (^) (= 0.70 MeV for a s = 0.4 and 4.2 MeV for a s = 0.8). The 
dashed lines relate to the same normalization condition (|3.18|) , but for \i = 5 GeV. Now, 
uo d(^) siny)£)( / u) ~ 0.1 MeV and the improper behavior of M n (m) has almost disappeared. 
However, if we use the normalization condition ( |3.23|) , we obtain lim m ^o^7r(^) = for 
arbitrary choice of /i (solid lines in Fig. 0). This particular scheme is therefore preferred, 
and we will discuss our results for the pion decay constant F n and the 7r° — > 77 decay width 
using this scheme. 

Finally, let us remark that from the low-momentum behavior of the SD solution (p one 
can read off the constituent quark mass m* as defined by Adler and Davis [0] . Reference J7| 
implies m* = — 1/<//(0). m* is a linear function of m and is depicted in Fig. |j. These results 
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are consistent with the constituent mass of Ref. [0] which is m* = 70 MeV for the massless 
case and the pure linear potential with o = (350 MeV) 2 . As a function of a, m* grows as 
square root (Fig. ||). 



V. THE PION DECAY CONSTANT F w 

So far we have been concerned with describing the hadronic structure which is determined 
in this context by the quark-quark interaction K. Therefore, for simplicity, we have so far 
omitted the external local operator L(x), which is used in the present approach to describe 
weak and radiative decays. However, rederiving the bilocal action W e ff in the presence of 
L(x) shows that we can consistently reinstate L(x) by the substitution 

M(x,y) -> M(x lV ) + L(x)6^(x -y) . (5.1) 



The matrix element for the leptonic decay of pseudoscalar mesons is then |jI7 | 



<l ± v l \W e f f \'K ± > = < l ± ui\^N c Ti(S(M + L)) 2 |7r ± > (5.2) 
= < Z ± i/ i |iiV c Tr(5jM5L)|7r ± > . (5.3) 

It is expressed through the axial-current matrix element which is conveniently parametrized 
by the pion decay constant F n . Evaluating (|5.3|) thus yields [[0|] for equal u and d quark 

masses 



In Ref. [IT] we obtained the first correct numerical results for the decay constants in the 



bilocal approach with the harmonic potential not only for the pion, but also for the kaon 
and their radial excitations. (Of course, (|5.4|) was generalized for different quark masses 



because of the kaon.) 

The variation of F w with o (for m = 7 MeV) is given in Fig. |^. grows monotonically 
with a roughly like a square root, Fig. [7]. 
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At n = 1 GeV and for fixed o = (350 MeV) 2 , our M v is fitted to the experimental value 
for m = 1.92 MeV when a s = 0.4, and for m = 2.80 MeV when a s = 0.8. Unfortunately, 
Fig. ^| shows that F n is then too small, being typically about 20 — 30 MeV . We remark 
that in the earlier treatments of the linear or funnel potential F n has been even smaller, 
e.g., F n = V2f w = 16 MeV in Ref. 0, and F n = (16 - 34) MeV in Ref. §. However, 
very close to the massless regime, at m = 0.22 MeV for a s = 0.4 and at m = 0.41 MeV 
for a s = 0.8, the correct experimental ratio M n /F n is obtained. This means that we can fit 
both M n and F v provided we rescale all dimensional quantities, including /x. Therefore, if 
we rescale by a factor of 7.9, we get M n = 140 MeV and F^ = 132 MeV at m = 1.71 MeV, 
a = (7.9 x 350 MeV) 2 at the renormalization scale fi = 7.9 GeV. Similarly, for a s = 0.8 
we again practically reproduce the experimental values if we increase the scale by a factor 
of 6: at n = 6.0 GeV, m = 2.47 MeV, and a = (6.0 x 350 MeV) 2 , we get M w = 140 MeV 
and F n = 132 MeV. Clearly, continually varying a s and o would give experimental values 
of M T and F n for continuum of different quark masses m but also different scales \i. 

VI. THE vr° -> 77 DECAY WIDTH 

If the external operator is taken to be L(x) = QA M (x)7 M , where Q = e diag(Q u , Qd, Q s ) = 
ediag(2/3, —1/3, —1/3), it will make possible the radiative processes as the n° — > 77 decay 
computed in [|17| for the harmonic interaction. We consider the present case of the funnel 
potential more realistic but still extremely oversimplified, so that we do not intend that 
our calculation should compete with the standard description via the Adler-Bell-Jackiw 
anomaly and PCAC, which yields an almost experimental decay width. We calculate this 
decay more as a further test of the quality of the funnel interquark interaction. We should 
however stress that our approach is in a way more ambitious than most of the other ap- 
proaches, including the standard anomaly calculation: these, namely, always contain the 
step when one actually parametrizes the unknown hadronic structure with the pion decay 
constant F n . On the contrary, the present calculation is not parametrizing but trying to 
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describe the pion structure and in this respect it is more microscopic. 

The 7T° — > 77 transition is caused by the cubic term from W e ff because it contains 
subterms with one meson bilocal M. and two photon fields A^. The transition matrix 
element is thus 

Xo 77 =< j(k,a)j(k',a')\iN c Tr [MSQASQfiS] |tt°(P) > , (6.1) 



where the symbol "Tr" also includes the integrations over coordinates. Equation (|6.1|) in 
fact corresponds to the 75 triangle graph except that the propagator lines emanate out of a 
pseudoscalar bilocal bound-state vertex and that these propagators are not free but dressed 
ones. Transforming to momentum space, 



T^o 77 = 2»iV c e 2 g " g % (fc, a)e v {k', o')I» v , (6.3) 

where 

^ = J ^- 4 tr[T(q ± \P)S(q - P)>fS(q + k')^S{q)] . (6.4) 

Inserting V and S, rearranging, and integrating over the parallel component qp and 
performing the spinor trace, 

7^ = Ae a ^ v — h , (6.5) 

® eI P = / 7Z-TaM<l±,k ± ,<p) — 2 ^- , (6.6) 

$tml fi = -lf J^M<1±, *4> <P){S[E(9±) + E ((<1 + k')±) + ^] 

x + L2(q±)} - 6[E(Q±) + E((q + k') ± ) - " Ho±)\} . ( 6 - 7 ) 
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Jp(q ± , k' ± , <p) = - f {q \ k ^ sin^(g ± ) cos^((g + k')±) 

+ i^k sin ^((g + k ') ± ) cos <p(q ± ) . (6.8) 

Since Ip is a function of M T and of only one four vector, (k' ± ), 

b = (k' ± )pC[E, L u L 2 , cp, , (6.9) 

where C is a dimensionless Lorentz scalar functional of E, L±, L 2 , and <p and a function of 
M„- and k' ± . We can thus extract it numerically by evaluating ( |6.9j ) in, say, the rest frame 
as the easiest one. Noting that e a/3flL/ P a (k' ± ) p = —e^^kak'^, summing l^o^) 2 over the 
polarizations a, a', and integrating over the phase space of the outgoing photons yields 

r(7r° -> 77) = a 2 M^\C\ 2 . (6.10) 

The dependence of Y on m (for fixed \i = 1 GeV and cr = (350 MeV) 2 ) is depicted in 
Fig. H The experimental value T exp = (7.7±0.5) eV is reached for a s = 0.4 at m — 11.7 Mel/ 
and for a s = 0.8 at m = 12.1 MeV. For such quark masses, M n is already too large, while F n 
is still too small. Since V falls rather quickly with m, in a good approximation proportional 
to m 3//2 , it is far too small in the range of m's for which we could fit both M w and F w using 
rescaling, as shown in the preceding section. 

On the other hand, for fi = 6.7 GeV, m = 25.6 MeV, and a s = 0.4 (or /i = 5.0 GeV, 
m = 23.3 MeV, and a s = 0.8) we can also fit F n and T to experiment using rescaling, but 
then M n becomes too large by a factor of 7 to 10. The variation of a with fixed m = 7 MeV 
does not yield better results either. This is reminiscent of the situation in the harmonic- 
oscillator case ]17| where it was impossible to fit M n , F„ and T simultaneously. In this 
case, however, we have not explored the whole parameter space and it is still possible that 
appropriate choice of a s , m, a, and \i would fit all three quantities simultaneously. 

The fact that T depends on the current quark mass as m 3 ^ 2 , is another manifestation of 



consistency with PCAC, as discussed in Sec. V of Ref. JT7 ] . 
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VII. CONCLUSION 



We have studied the pion as the quarkonium bound by the funnel potential. We did 
it using the covariant generalization of the instantaneous-potential model proposed in the 



framework of the effective bilocal Lagrangian [13-T6|]. Provided that the covariant BSE 
solutions thus obtained describe bound states sufficiently well, they can be useful for un- 
derstanding what happens in the experiments probing the quark substructure. In these 
experiments both nonperturbative bound-state effects and relativistic recoil effects can be 
important, as in the program of CEBAF, for example. 

Our work has extended the line of research which employs instantaneous interquark 
potentials (e.g., Refs. ||-§,[10],[TT| and [031 ) by using the boost-invariant potential ansatz, 
by introducing and analyzing the effects of the finite quark mass, and by improving and 
generalizing the subtraction procedure for UV divergences caused by the Coulomb part of 
the interaction. The initial choice for the renormalization prescription, essentially following 
from the standard prescription where infinities are subtracted at a spacelike point |26|, has 
been shown to preclude the correct chiral limit. However, by appropriately modifying the 
renormalization prescription, we have been able to recover the correct chiral limit. With this 
(naturally preferred) scheme, our pion exhibits the qualitatively correct (pseudo)-Goldstone 
behavior. We have also shown how the MP breaking occurs in this context and explored 
the dependence of the results on the subtraction point //. In particular, we have found that, 
as fi grows, the results in the first renormalization scheme gradually approach the results in 
the second, preferred scheme. 

As far as the quantitative results are concerned, the description of the pion is still not 
satisfactory, as the usage of the funnel potential has not (yet) resulted in the improvement 
of M n , F n , and r(7r° — > 77) with respect to the values we obtained with the harmonic 



interaction [17 |. 

These results are not definitive as we have not yet systematically explored the parameter 
space. This may be one task for the future work, but further search for an improved form 
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of the potential certainly remains as the other, and even more important task. 
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APPENDIX: 

The SD equation ( |3.26| ) and the BS equation (|4.4|) are nonlinear integral equations with 
singular kernels. In order to solve these equations numerically, we discretize them, i.e., we 
approximate integral equations with a finite set of coupled nonlinear equations. This is 
realized by discretization of the momentum variable, 

h = k max C-^) n (i = l,...,N + l) . (Al) 

kmax is the "numerical" cutoff; it should be chosen large enough to eliminate the effects 
of the boundary condition at "infinity" . N + 1 is the number of points. For n = 1 Eq. ( |A1| ) 
corresponds to an equidistant mesh, i.e., uniformly distributed points, while for n > 1 the 
points are denser near k = 0. A finite set of variables is defined, e.g., for the SDE we define 
(pi = <f(ki), i = 1, N + 1. The values (p\ and ^n+i are fixed by the boundary conditions 
f(0) = 7r/2 and lim fc ^ +00 ip(k) = (see, e.g., |J). Equation ( |3.26|) must be satisfied for 
every k i: i = 2, N; thus we have obtained a system of N — 1 nonlinear coupled equations 
for if i, ...,</jjv_i. 

A problem arises from the singular nature of the integration kernel. The integrals have 
to be calculated as principal values. So, we define the continuous function <f(k) as a cu- 
bic spline of <pi, <pn+i and perform integration using Gaussian quadrature, adapted for 
principal-value calculation. The system of N — 1 coupled nonlinear equations is solved us- 
ing a modification of Brent's methods. A similar procedure was applied to the BSE ( |4.4[ ), 
treating M n as an additional variable. 

22 



For k max larger than a few MeV, the solutions has been found to be independent of 
k m ax- For an equidistant mesh, n — 1, the number of points have to be N > 100 to obtain 
a solution independent of N. For a mesh with n — 3, this is TV 25. The meshes with 
larger iV have an inadequate distribution of points in respect to the solution (p(k). This 
optimization is important because the computer time consumed behaves as oc iV 2 . 
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FIGURES 

FIG. 1. Solutions <p of the SDE, Eq. Q), and L x and L 2 of the BSE, Eqs. flOP , are plotted 
versus q for m = 7 MeF and a = (350 MeF) 2 . Solid lines pertain to the case a s = 0.4, and dashed 
lines to a s = 0.8. Starting values of the solutions are ip(0) = 7r/2, Li(0) = 1, and £2(0) < 0. 

FIG. 2. Pion mass vs. quark mass m for three different values of a s (0.0,0.4,0.8) and 
a = (350 MeV) 2 . The dotted lines essentially correspond to the normalization condition ( p. 18 ) 



and fi = 1 GeV . The dashed lines relate to the same scheme, but for /j, = 5 GeV. The solid lines 
are obtained by using the normalization condition ( 3.23| ), slightly modified, as explained in Sec. ID. 
The renormalization point is fj, = 1 GeV. 

FIG. 3. Pion mass M w vs. string tension a. As in Fig. ||, the solid lines relate to the normal- 
ization condition ( |3.23| ). 

FIG. 4. Constituent mass m* vs. quark mass m. The curves are as in Fig. |2[ 

FIG. 5. Constituent mass m* vs. string tension a. The curves are as in Fig. ^. 

FIG. 6. Pion decay constant F n vs. quark mass m. The curves are as in Fig. 

FIG. 7. Pion decay constant F w vs. string tension a. The curves are as in Fig. ||. 

FIG. 8. Decay width for 7r° — > 77 vs. quark mass m. The curves are as in Fig. ^. 

FIG. 9. Decay width for tt° — > 77 vs. string tension a. The curves are as in Fig. 0. 
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